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NOTE SUR LA LOI DE RÉSISTANCE DE L'AIR 


ET SUR LA CORRECTION D'ÉLASTICITÉ 
PROPOSÉES PAR M. LE CAPITAINE DARRIEUS. 


Sur la Semande qu’a bien voulu me faire M. le général Mochot, j'ai examiné 
au point de vue théorique et soumis au contrôle de l’expérience les idées émises 
par M. le capitaine Darrieus dans ses Notes datées du 20 février et du 5 mars 
1918 «sur une nouvelle cause de variation de portée avec la température ». 

Les ésultats que j’ai obtenus confirment entièrement les prévisions contenues 
dans ce travail que je considère comme étant de tout premier ordre, tant par 
l'importance des indications nouvelles qu’il apporte que par la pénétration 
dont l’auteur a fait preuve en établissant par des considérations simples, mais 
. discutables, de théorie cinétique, une forme de la loi de résistance de lair qui 
se trouve être exacte aux vitesses des projectiles d'artillerie, mais qui cesse de 
l'être aux faibles vitesses lorsque les effets de viscosité ne sont plus négligeables 
par rapport aux effets d’inertie et d’élasticité du milieu. 

Je considère, en particulier, comme tout à fait remarquable la manière dont 
le capitaine Darrieus déduit, par application d’un théorème de similitude balis- 
_ tique, les conséquences pratiques de sa loi de résistance, et montre comment 
les coefficients de correction déjà connus perméttent d’exprimer la correction 
nouvelle. 


I. Similitude aérodynamique et loi de résistance. 


Voyons tout d’abord comment tl est possible de déduire des qe générales 


de l'aérodynamique la forme de loi proposée par M. Darrieus, sous la condition que ` 


les termes de viscosité puissent être négligés. 

._ Le résultat énoncé par M. Darrieus est le suivant :`si on désigne par v la 
vitesse du mobile, par a la vitesse du son dans le fluide de densité À, où le 

mobile se meut et par æ une dimension linéaire de ce mobile, la résistance F 

à laquelle il est soumis est de la forme : 


(1) . FA (©) 


où la fonction Ÿ ne dépend que de la forme du projectile. © ~ 

Notre démonstration sera simplifiée si nous mettons la relation (1) sous une 
autre forme également utile pour la vérification expérimentale dont il sera 
question plus loin. 

Remplaçons dans (1) la densité À, par sa valeur tirée de la formule de La- 
place : 
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ôli 
où p, représente la pression du fluide dans les régions où la perturbation due 
au passage du projectile n’a pas encore ps On obtient la forme : 


6) Foei) 


| NC) 


C’est la relation (3) que nous nous proposons d’établir pour le cas du régime 
permanent de mouvement de translation d’un projectile dans un fluide com- 
pressible non visqueux. 

L'absence de viscosité implique l'existence à tout instant, en tout point du 
fluide, dans la région troublée par le passage du mobile, d’une pression iso- 
trope que nous représentons par p au point de coordonnées x, y, z, et à l’ins- 
tant t. 

Sur chaque élément de surface ds du mobile, le fluide exerce, en vertu de 
la même hypothèse, une action élémentaire normale égale à p ds. 

L'ensemble de ces actions se réduit, dans le cas général, au système d’une 
force et d’un couple résultant. La composante de la force dans la direction 
d’un axe, celui des x par exemple, a pour valeur : 


(4) F —[r cos B ds, 


B étant l'angle que fait, avec laxe des x, la direction normale à l'élément de, 
et l’ intégrale étant nai à toute la surface du projectile. 

La manière dont la pression p varie autour du projectile en mouvement 
est déterminée, par rapport à celui-ci, dans le cas du régime permanent, par 
les équations habituelles de l’aérodynamique où Vz, Vy, v, représentent les 
composantes de la vitesse au point x, y, z, êt A la densité : 


1 dp dv, dv, dv, 
dd 13e FH Pr 
1 dp dy dv, 
5) A à Va So Yd “ia "s Ja 
1 dp dv, x dv, 
er Va a + 19 PY 19 


auxquelles il faut joindre la condition de continuité : 


(6) | ae a "ea +20 


ainsi que les conditions aux limites : 
1° Vitesse normale nulle sur la surface du mobile ; 
2° Vitesse égale à — v à distance infinie, 


20, 
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Enfin la pression p et la densité p sont liées par la loi de compressibilité 
que nous pouvons, en raison de la rapidité du mouvement, supposer adiaba- 
tique, et qui s’ exprime par les conditions suivantes : 

3° En dehors des surfaces de discontinuité, comme celle de l’onde de choc, 
par la condition : 


à d d\ p 
(7) | (LÉ ete vieu 


4° Sur la surface de l’onde de choc, en un point où la normale à celle-ci 
fait un angle 8 avec la direction de la vitesse v, par la relation d'Hugoniot : 


NT pe Re 
(8) ae nn Re v? sin? 0. 


I s’agit de montrer que toutes ces conditions satisfont à un théorème de 
similitude que nous établirons de la manière suivante : 
Introduisons des variables réduites définies par les relations : 


E: atoy Oz 
E=- n=% € 
do — = ds 
v, v v, 
v 
Ue 


L'expression de la composante F devient : 
(10) F = a°p, | w cos do: 
La distribution de la pression réduite & est déterminée par des équations 


“et des conditions faciles à obtenir par substitution, et qui se présentent sous 
la forme suivante : 


4 ds du, du. 
SET. RE 
FES 1 dw du 
(5) Sn a, Su Sr + u z na 
1 ds du, 
7 K ny + tu e 
: ; d (ðu) , d(du,) | d (du.) dd 
(6) ET UE LEE FS 


1° Valeur nulle de la composante normale de u sur la surface. 
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2° Valeur de a égale à — U à l'infini. } 
3° En dehors des surfaces de discontinuité : 
W (at onten)? 
4° Sur la surface de londe de choc, on a la condition réduite d’Hugoniot : 
(8) (y+ 1)w — (y — 1) = 2yU’ sin? 0. 


Ces équations ne font intervenir que les variables réduites et sont les seules qui 
doivent être satisfaites par la distribution du mouvement en régime permanent. Si 
donc F ai distribution de régime permanent est possible dans un gaz autour d'un pro- 
jectile en mouvement par une valeur donnée de U, une distribution homologue de à 
et ò est possible autour d’un projectile semblable se mouvant dans un gaz de même y 
avec la même valeur de U. 

Toutes ces distributions homologues correspondent à une même solution des 
équations réduites, caractérisée, en particulier, par une distribution déter- 
minée de w en fonction de Ë, », €, U et y, la forme de cette fonction dépen- 
dant uniquement de la forme du projectile et de son orientation par rapport 
à la direction de la vitesse réduite U. 

L'intégrale : 


[a c0s Bd = pu) 


est donc elle-même une fonction de U dont la forme dépend de y, des para- 
mètres caractéristiques de la forme du projectile et de son orientation par 
rapport au mouvement. 

En introduisant ce résultat dans l'expression (10) de la résistance F et en 


remplaçant U par s, on obtient le résultat qu'il s'agissait d'établir : 
F = a°p,Ẹ £) 


. Les trois composantes de la force exercée par le fluide ont des expressions 
de cette forme; on verrait aisément de la même manière que le moment du 
couple résultant a des composantes de la forme : 


(11) -ap ð (2). 


On peut encore énoncer les résultats précédents sous une autre forme qui 
est peut-être plus saisissante. 

La distribution du mouvement dans un gaz autour d’un projectile) ne dépend 
que de la forme de celui-ci, du rapport de sa vitesse à la vitesse du son et du rapport y 


des chaleurs spécifiques, 


+ 


—— 


(1) La similitude s'étend, naturellement, à la forme et à la position de l'onde balistique 
autour du projectile, lorsque la vitesse de celui-ci est supérieure à la vitesse du son. 
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Si nous étions partis des équations complètes de l’aérodynamique en y fai- 
sant figurer les termes de viscosité, nous n'aurions pas trouvé la même simi- 
litude.:J’ai rappelé plus haut, et une évaluation facile permet de le démontrer, 
que ces termes sont négligeables dans les conditions habituelles des tirs de 
l'artillerie. 

On peut comprendre le rôle important joué dans ces phénomènes par le 
rapport y des chaleurs spécifiques en faisant intervenir les mouvements de 
rotation des molécules ou leurs mouvements internes dans le raisonnement 
cinétique sommaire, mais intuitif, que fait M. Darrieus dans sa première Note 
et qui l’a certainement mieux servi comme instrument de recherche que n’aurait 
pu le faire une démonstration plus précise. | 

Bien qu’il ne s’agisse pas exactement d’un régime permanent analogue à 
ceux qui ont été envisagés dans la démonstration précédente, il est intéressant 
de vérifier le résultat de M. Darrieus dans un cas particulier pour lequel on 
peut conduire le calcul jusqu’au bout. 

Jl s’agit du cas d’un plan indéfini se mouvant dans la direction de sa nor- 
male, avec une vitesse supérieure à celle de la propagation du son. Il se produit 
de chaque côté du plan une onde d’Hugoniot, séparée du plan par une région 
comprimée à lavant et raréfiée à l'arrière, et dont la distance au plan augmente 
proportionnellement au temps. 

Il est facile, étant données la vitesse v du plan et les conditions du milieu, 
de calculer les pressions avant et arrière entre le plan et les deux ondes de 
choc, et d’en déduire la résistance au mouvement, proportionnelle à la diffé- 
rence de ces deux pressions. Le résultat est bien de la forme indiquée par la 
théorie avec la valeur suivante pour la fonction @ (U) : 


ERNETEN 


d’où pour la fonction y de la formule (1) la forme simple `: 


Ce résultat n’est valable que pour les valeurs de U supérieures à l'unité et 
donne pour représenter la fonction f(v), au-dessus de la vitesse du son, une 
courbe décroissante d’allure conforme à celle de la courbe de Gâvre, et tendant 
“Eur, 

CE 
Ceci correspondrait, pour le plan se mouvant normalement à sa direction, à 
un indice de forme voi in de 5 aux grandes vitesses, résultat qui semble par-- 
faitement acceptable a priori. > 


aux très grandes vitesses vers une asymptote horizontale d’ordonnée 
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II. Similitude balistique et corrections de tir. 


Si nous passons maintenant à l'étude du mouvement d’un projectlie sous 
l’action combinée de la pesanteur et d’une résistance conforme à la loi démon- 
trée plus haut, nous allons établir un théorème de similitude balistique, tout 
à fait indépendant de la similitude aérodynamique envisagée jusqu'ici et qui 
va nous permettre de retrouver, sous une forme un peu différente, les résultats 
de M. Darrieus sur la manière dont il convient de faire les corrections de tir 
pour tenir compte de la loi exacte suivie par la résistance de Pair. 

Si m est la masse du projectile, la retardation due à la résistance a pour 
valeur : 


(12) r=n—ho() 


c étant un coefficient analogue au coefficient balistique ordinaire. 

Les équations du mouvement du centre de gravité sont, en désignant par 1 
l'inclinaison de la tangente à la trajectoire sur l'horizon, et en supposant que 
la résistance n’a pas de composante transversale : 


(13) : 
: = RQ (U) cos r --g 


Envisageons d’abord le cas d’une atmosphère homogène pour laquelle p, 
et a sont constants, et introduisons des variables réduites définies par : 


| =" pi 

Fear so pet 

(14) 3 \ ET RES =; 
| K Po 


(15) 


Il suffit, pour assurer la similitude des trajectoires, que K soit constant, 
c’est-à-dire que le coefficient balistique varie en raison inverse A la pression ambiante 
et que la vitesse initiale varie proportionnellement à la vitesse du son. 

Dans ces conditions, la portée, qui se comporte comme les coordonnées d'espace, 
varie proportionnellement au carré de la vitesse du son. 
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La portée réduite : 
RE rs 


aÊ 


ne dépend que des deux éléments U, et K, pour un angle donné de tir, On a 
donc pour les corrections : i 


dr dU, dK 
en désignant par U, la vitesse initiale réduite : 
CRIN 
u= 


si V, est la vitesse initiale elle-même. 
En remplaçant dans (16) les variables réduites par leurs expressions, il vient : 


3 ðX dV, de  dp, da 
Un, Sr -r)+e 


Comme la vitesse du son ne dépend que de la température, si le gaz reste le 
même, et varie proportionnellement à la racine carrée de la température absolue 
© on peut écrire : 

dX dV, dc dp, C. \ 26, 
G Toena O E a 
où les corrections sont exprimées en fonction des indications lues directement sur le 
baromètre et le thermomètre. 

Si la composition du gaz change, y restant le même, on peut conserver comme 
variables la pression et la densité, et on retrouve exactement la forme donnée 
par M. Darrieus : 


Ce mode de correction, pour tenir compte de la variation du degré hygromé- 
trique, est cependant sujet à caution parce qu’il suppose constante la valeur 
du rapport y des chaleurs spécifiques et que ce rapport n’est pas le même pour 
lair et pour la vapeur d’eau. Il sera:très probablement suffisant, au point de 
pratique, de s’en tenir à la forme (18) qui est la plus simple et la plus immé- 

iate. 

Je rappelle que le coefficient de correction C,, lorsque la variation relative de 
vitesse le long de la trajectoire est faible, a pour valeur : 


v 


La fonction f est la fonction ay habituellement considérée et qui se con- 


fond, à un facteur constant près, avec la fonction ẹ de la formule (1). 
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III. Généralisation de la similitude balistique. 


Je me propose de montrer qu’il est possible d'étendre le théorème de similitude 
balistique au cas où la température © et la pression p varient avec l'altitude y et de 
trouver la condition nécessaire et suffisante à laquelle doivent satisfaire les 
variations simultanées de p et de © pour que le théorème subsiste et que la 
formule de correction (18) reste applicable, exactement sous la même forme où 
poet ©, désigneront maintenant la pression et la température au sol. 

Si A est la densité de l’atmosphère au point d'altitude y, on a, en désignant 
par R la constante des gaz parfaits pour l'unité de masse do cette atmosphère : 


(20) p= ROA. 


a désignera toujours la vitesse du son, fonction de l'altitude y et donnée par la 
formule : 


(21) a =4/ y = VRO. 
La statique des gaz donne en outre la relation : 


(22) dp =— Ag dy. 


v = Vv? + v? désignant toujours la vitesse du projectile, les équations du 
mouvement sont : 


(23) 


Introduisons des variables réduites définies par les relations suivantes : 


dx — adË dy = a°dn dt =a dû 
dé vy d 


v Ur + oyati 
i ss AT ee 


d9 Jid 


(24) 


a étant fonction de y, l équation différentielle dy = adn définit n en fonction 
de y, à une constante arbitraire près que nous supposerons choisie de manière 
que » soit nul en même temps que y à l'origine de la trajectoire, où l’état du 
milieu est caractérisé par po et Oo. | 
L’équation différentielle (22) devient, en tenant compte de (20) et (21) : 
dp 
ne di: 
p gy an 
d’où : 
—gyn. 
(25) D = pol I 
h Zags 1 700. Í Se 
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La variable y est donc liée de manière très simple au mepet 4 - dont elle 


1 


est une fonction représentée por : 
~ 4 p 
(26) = — — log =. 


Les équations du mouvement (23), exprimées au moyen des variables réduites 
. deviennent, en tenant compte du fait que a est fonction de y : 


zU us aiee TIA i COS T 
(27) oy c 


a Po, 971 . 
S GT re @(U)sin r—g. 


Pour que ces équations se conservent quand on passe d’une trajectoire à 
une autre correspondant à de nouvelles conditions atmosphériques, il ēst 
nécessaire : x 


1° Que 2° reste le même, FE ER que le coefficient balistique varie 


propor onelaiti à la pression au sol; 
2° Que la fonction de. n qui représente — ga reste la Lo qu’elle soit une 


fonction permanente, indépendante de l’état fe l atmosphère. Z 
Nous pouvons mettre cette condition sous une forme plus simple et plus 
intéressante. 
En tenant compte de la relation : 


dy — «°dn — RO dn 
la condition s'écrit, si la fonction f désigne une fonction permanente : 


2 = f(n) d 


d’où, en tenant compte de (26) et en représentant par F une autre fonction 
permanente : 


(28) a —F(E) 


Ce résultat est remarquable. La pression et la température peuvent varier 
de manière quelconque avec l'altitude, sans que le théorème de similitude 
cesse d’être exact, à condition que les variations simultanées de la pression et 
de la tempétature, quand on pen d’une trajectoire à une autre, soient telles 


que à 2 reste une même fonction de 2 i quelle que soit d’ailleurs cette fonction. 


Ce” résultat général comprend, comme cas particulier, le cas d’une atmo- 
sphère à température uniforme où la pression varie avec l'altitude suivant la 
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£ 
loi de Laplace ene " puisque, de quelque manière que varie cette température 


Rae © ; ; 

uniforme d’une trajectoire à l’autre, o est bien une fonction permanente 
Là s x TS à g i 

égale toujours à l'unité. 


Un autre cas particulier, le plus intéressant de beaucoup, parce qu'il corres- 
pond en moyenne exactement à la réalité, est celui où la variation de tempéra- 
ture avec l'altitude correspond à la loi de détente adiabatique, ainsi que le 
signale très justement M. Darrieus à la fin de sa Note du 20 février. 

En effet, on a, dans ce cas, suivant la loi connue : 


29) AE 


Cette relation est bien Ja forme (28); la fonction permanente F étant une 


puissance dont l’exposant Ee dépend uniquement de la nature des gaz 
qui composent latmosphère. 

Enfin, pour que la similitude entre les trajectoires existe, il ne suffit pas que 
les équations différentielles du mouvement exprimées au moyen des variables 
réduites conservent leur forme, il faut encore que les constantes d'intégration 
réduites soient les mêmes, par conséquent langle de tir et la vitesse initiale 
réduite U,, c’est-à-dire que la vitesse initiale elle-même V, doit varier propor- 
tionnellement à la racine carrée de la température au sol VOo. 

On voit facilement que les conditions de similitude qui viennent d’être recon- 


nues nécessaires sont en même temps suffisantes. La variation de portée s'obtient 
immédiatement en tenant compte de (26) et (28) : 


(30) PS f de= f èdE— yR f QdE— à f © (n)dé, 


® étant une fonction permanente reliée de manière très simple à la fonction F 
par une exponentielle. 


La portée varie donc comme la texapérature au sol quand les conditions de 
similitude sont remplies. | 


Nous sommes ainsi conduits au remarquable énoncé suivant : 


Si l'état de l'atmosphère change d'un tir à l'autre, de manière que le rapport de la 
température (absolue) en un point quelconque à la température au sol reste une même 
fonction quelconque du rapport de la pression en ce point à la pression au sol, la portée 
varie comme la température au sol lorsque, l'angle de tir restant le même, le coefficient 
balistique varie comme la pression au sol, et la vitesse initiale comme la racine carrée 
de la température au sol. | 

Les formules de correcuon de tir proposées par M. Darrieus restent exactes sous la 
seule condition très générale, que l’état de l'atmosphère change de la manière indiquée 
et en particulier en suivant toujours la loi adiabatique, conforme, en moyenne, 
aux faits d'observation météorologique. 
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Les formules de correction ont ceci de remarquable qu'elles ne font intervenir 
que des coefficients déjà connus et des observations faites au sol. 

Il y a lieu, à mon avis, de les adopter dès maintenant, tant à cause de leur exactitude - 
théorique que des vérifications expérimentales auxquelles j'ai pu soumettre la nouvelle 
loi de résistance de lair: : 

Avant d'indiquer le principe et de donner les premiers résultats de ces 
recherches, il peut être intéressant de voir quelle forme particulière prennent, 
dans le cas de la distribution adiabatique, les fonctions qui figurent dans la 
théorie générale de la similitude. E 

La relation caractéristique : 


jointe à (25) donne: 
d'où : 


et par intégration : 


ANR iig 
O0 =0,— eT y. 

La température varie en fonction linéaire de l'altitude, mais ce n’est qu'un 
cas particulier et cette relation n’est nullement nécessaire, ainsi qu'on l’a vu, 
pour que la loi de similitude et la formule de correction qu’elle entraîne soient 
exactes. : We 

La vitesse du son varie en fonction de l'altitude suivant la relation : 

a = yRO, — (y — gg. 


En vertu de la relation (30) la portée est donnée par : 
X—yR@, fet = În qé. 


La relation entre Ë et n permettant d'effectuer cette intégration dépend de 
la forme admise pour (U), et s’obtiendra par intégration des équations (27). 
qui deviennent, en supposant les dérivées prises par rapport à 7 : 


UE ey he 7 (VE Fn) co 


— 1 ʻ EEE N 
j reg : y2 he YM o (Vg Fa?) sin 7—g. 
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On voit facilement, èn tenant compte des relations introduites entre les 
variables, que la condition (28) peut encore s’écrire : 


G n) 


F, étant une fonction quelconque. Cette forme a l'avantage de ne faire inter- 


- venir que les données des sondages de”température. L’énoncé du théorème 


général de similitude devient : 


Si l’état de l'atmosphère change d’un tir à lautre de manière que le rapport de la 
température (absolue) en tout point à la température au sol reste une même fonction quel- 
conque du quotient de l'altitude du point par la température au sol, la portée varie comme 
la température au sol lorsque, langle de tir restant le méme, le coefficient balistique 
varie comme la pression au sol et la vitesse initiale comme la racine carrée de la tempé- 
rature au sol. 


Les formules de correction (18) et ( 19) restent applicables dans ces condi- 
tions. 


IV. Vérification expérimentale. 


Les recherches que j'ai entreprises en collaboration avec MM. Vaillant et 
Saphores sur l’utilisation des courants d’air à grande vitesse pour la solution 
des problèmes de balistique, trouvent dans la question actuelle une application 
immédiate et simple. 

Supposons en effet que la veine gazeuse soit obtenue par détente dans 
l'atmosphère (de pression pọ) de l’air comprimé contenu dans un réservoir sous 
la pression p.. 

Les formules classiques montrent que le rapport de la vitesse v acquise par 
lair détendu à la vitesse a du son dans les conditions où se trouve cet air ne 
dépend que du rapport z par la relation : 


0 


ga 
yl B) y ah: 
a AA BNA 
Si la loi proposée est exacte, la poussée exercée par ce courant d'air sur un 
obstacle de forme quelconque est donnée par : 


F = kpç® (2) ; 


où k est une constante qui dépend uniquement de l'obstacle. 
Po étant constant, la force observée ne doit dépendre que de la pression p, quelle 
que puisse être la température de lair comprimé dans le réservoir sous cette pression. 
Or il nous est facile de faire varier cette température. Le réservoir étant 
rempli à la température ambiante, sous diverses pressions de remplissage P 
que nous avons pu faire varier de 2 à 6 kg, la température de l'air qui reste dans 
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ce réservoir au moment où sa Re est tombée pendant l’évacuation à une 

même valeur p, est d'autant plus basse que la pression de remplissage a été plus 
élevée. Par exemple, pour une pression p, de 2 kg, la température est ambiante, 
20° par exemple pour P = 2 kg et égale à — 60° environ pour P = 6 kg. 
Nous avons ainsi une marge de température de 80°, que nous pourrons porter 
au delà de 300° quand l'installation en cours sera achevée. 

Le mode expérimental est par conséquent très simple : nous enregistrons les 
efforts exercés sur l'obstacle pour une même valeur de la pression dans le réservoir et 
diverses valeurs de la pression de remplissage. Ces efforts doivent rester les mêmes 
bien que la vitesse et la température de la veine varient dans de rh limites. 

Les vitesses ont varié de 100 à 400 m par seconde. 

Les nombres ci-dessous sont relatifs à une tige cylindrique de 6 mm de 
diamètre. 

D’autres enregistrements ont été pris pour divers obstacles et M. Vaillant en 
fait actuellement les lectures. 

En face de chaque pression p, inscrite dans la première colonne du tableau 
se trouvent inscrits en ligne horizontale les efforts mesurés F relatifs à des pres- 


sions de remplissage de 6, 4 et 2 kg respectivement. La loi prévoit la constance 


des nombres dans chaque ligne horizontale. 


E-z 

Pr EEG F2 P=2 

1,06 9,44 0,43 0,42 

1,12 0,83 0,82 0,83 

KISS :. _1,32 1,29 1,28 

1,37 2,01 2,04 2,03 

1,44 2,22 2,22 SEED OR 

1,51 PT E9 53 Shen 2,53 
k 1,58 2,89 2,85 2,85 

1,65 3,22 3,18 | 3,29 

1,72 3,58 3,50 3,13 

1,87 4,09 3,95 4,03 

1,95 4,51 4,49 - 4,48 

tt 5,87 5,91 

2,19 6,51 6,46 

2,28 7,05 7,00 

2,37 7,53 7,53 


La concordance est excellente et les divergences sont de l’ordre des erreurs 
expérimentales. 
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